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To evaluate the feasibility of the recently proposed detection scheme of low energy neutrinos 
released from the Sun and supernovae called CLEAN, Cryogenic Low Energy Astrophysics with 
Noble Gases, which relies on the transparency of noble-gas cryogenic liquids to VUV radiation 
produced by neutrinos, we analyze theoretically VUV emission, self-absorption, and scattering of 
liquid helium and neon, primary candidates for CLEAN. Owing to strong repulsion of noble-gas 
atoms in the ground states at the equilibrium distance of the relevant excited state, the emission 
spectrum is substantially shifted from the absorption spectrum, and in principle the absorption 
is expected very small, allowing building large detectors. Our analysis, however, shows that the 
self-absorption and Rayleigh scattering are comparable to the size of the proposed detector. 

Our theoretical emission spectra are found in agreement with experimental observations although 
some deviation exists due to binary-interaction approximation, and our ab initio Rayleigh scattering 
lengths are found in agreement with other calculations based on the extrapolation of experimental 
refraction indices. The absorption process can result in either re-emission, which conserves the 
number of photons but delays their escape from the liquid, or in non-radiative quenching. 

PACS numbers: 32.30.Jc, 32.70.Fw, 29.40 Mc, 31.15.Ar, 31.70.-f 



Recently, McKinsey and Doyle pj has proposed the 
detection of low energy neutrinos released from the Sun 
and supernovae with the scheme called CLEAN, Cryo- 
genic Low Energy Astrophysics with Noble Gases. The 
central element of a CLEAN detector is a 6-meter diam- 
eter Dewar filled with liquid neon which is functioning 
both as a passive shielding medium and an active self- 
shielding detector. The advantage of liquid neon is high 
scintillation yield and low radiactive background. Other 
noble-gas liquids can be used as well, of particular in- 
terest is liquid helium. The self-absorption and scatter- 
ing of UV light produced by these liquids are critical for 
determination of the size of the detector, the interpre- 
tation of neutrino events, and exclusion of background 
events. The emission of liquid helium and neon in UV 
has already been studied E EI and scattering lengths 
of noble-gas liquids at scintillation wavelengths has been 
estimated from extrapolations of refraction indices 
the self-absorption is yet to be measured. Because there 
are many factors such as spectral reshaping and the ef- 
fects of impurities (for example, 0.2% of N2 impurities in 
liquid helium reduces the emission intensity by 20-40% in 
the region 700-850A Q) that influence the interpretation 
of the experimental results, the theoretical calculations 
of emission, self-absorption, and scattering, which con- 
stitute the subject of this letter, are important. Spectral 
properties of noble-gas liquids are not well understood 
mostly due to interdisciplinary nature (atomic, molec- 
ular, condensed-matter physics) of the phenomenon as 
well as the scarcity of the data, so this paper also has 
the aim of attracting the attention of theorists of diverse 
expertise and motivating further UV spectroscopic mea- 
surements. 

The scintillation of liquid helium due to energetic a- 
particle stopped in the liquid was first studied by Moss 



and Hereford [|| and by Kane et al. 0. Moss and Here- 
ford observed significant decrease of the intensity be- 
low A-point, and attributed this to the change of the 
oscillator strength of 2P-1S He transition due to the dif- 
ferences in spatial distribution of nearest neighbors in 
the superfluid and the normal fluid. The emission spec- 
trum of liquid helium was also studied experimentally 
by Stockton et al. 3] who obtained and compared UV 
spectra of electron-bombarded superfluid helium at 1.4 
and 2.1 K with those of high-pressure helium discharge 
lamp in the range of 1000-600 A. The spectra were very 
similar, but the liquid spectrum had an earlier cut-off at 
about 600 A due to stronger absorption. Keto et al. Q 
and Hill et al. 9] showed that the ionization of the he- 
lium results in formation of He^A 1 ^^) and He^a 3 !]^) 
dimers in the liquid which decaying produce UV scin- 
tillation. A number of publications [loL fill IT3 . IT^ ad- 
dressed an unusual phenomenon in liquid helium such as 
the formation of cavities around free electrons and ex- 
cited atoms and molecules. The dynamics of the forma- 
tions of cavities is fast, 5 x 10 -11 s 0, so that excited 
molecules do not decay significantly during the time of 
cavity formation in which case the liquid effects on the 
spectral properties are greatly reduced. 

The luminescent spectrum of liquid and solid neon was 
observed by Packard et al. £j in the range 700-900 A. It 
is interesting to note that the solid has a narrow peak 
at 743 A , but the liquid does not; similar peaks in solid 
and liquid exist at 774 A with widths about 30 A. The 
luminescence of neon gas was also briefly investigated: 
at low pressure P < 40 Torr, a single narrow line was 
detected near 743 A. At higher pressures, the intensity 
of the 743 A line decreased but the second line at 800 A 
appeared. At low temperature, 25 K, the 800 A line does 
not persist. 
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Absorption and emission spectra of liquids can be ob- 
tained approximately from the spectra of molecules. For 
example, in Ref. 3] it was found that lines in the 600A 
band of He2 appear in the liquid emission spectrum and 
can be well predicted from potential curves of the pair in- 
teraction using Franck-Condon principle. We would like 
to apply this principle to obtain theoretical UV emis- 
sion in helium and neon liquids (molecules) near 16 eV. 
The main peak of helium scintillation according to Ref. 3] 
appears around 15.5 eV (800A), and neon liquid scintil- 
lation was observed at 16 eV |4|- For calculations of the 
wavelength of the maximum of emission it is enough to 
know the equilibrium interatomic distance {R e q) of the 
lowest singlet excited molecular state, the energy of this 
state, and the interatomic energy of two ground state 
atoms at the distance R eq . Detailed spectral distribu- 
tion can be obtained from the overlap of wave functions 
of nuclear motion, according to the classical work by Con- 
don j 1 ll ] . Absorption transitions start from ground state 
atoms which are separated by a relatively large distance. 
Thus the small energy shift due to environment can be 
neglected and it is a reasonable approximation to con- 
sider only the strongest allowed atomic transitions which 
have transition energies close to the emission energy of 
about 16 eV. 

The helium interatomic or molecular potentials are 
well-known 

El El IH El • The potential curves of the 
lowest excited states of Ne2 are given in Ref. Ell- The 
ground state of neon diatomic molecule is investigated 
in [2I El I2I, . Analytical expansions of the ground- 
state interatomic potentials of noble-gas atoms that re- 
produce well the available vibrational spectra are given 
in [2i|. Most studies do not consider small interatomic 
distances in detail which are important for our applica- 
tion, so we supplement them with our own calculations 
described below. 

Our ab initio ground-state interatomic potentials at 
small distances are calculated using Pauli's principle. 
When two closed-shell atoms have wavefunctions over- 
lapping, electrons are promoted to excited atomic states 
increasing the energy of each atom. In helium, the repop- 
ulation increases the energy by 22 eV times overlap and 
in neon 16 eV times overlap. Wavefunction overlaps are 
obtained using Dirac-Hartree-Fock potentials, Is 2 in the 
case of helium and ls 2 2s 2 2p 6 in the case of neon. Over- 
laps of all electrons of one atom with those of the other 
are taken into account. In the range of interatomic en- 
ergies of interest for self-absorption other effects such as 
electrostatic direct and exchange interactions are weaker 
because of screening. We have checked the magnitude of 
electrostatic effects in case of helium, and found that the 
dominant contribution comes from the Pauli principle. 

In FigQ we show the comparison of the ground-state 
helium interatomic potentials obtained from the ove rlap 
theory, the extrapolation using the expansion of Ref. [24| , 
and ab initio calculations with exponentially correlated 
Gaussian (ECG) functions E3- The expansion of [24j 
agrees well with calculations of E3 up to 10 eV indicat- 
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FIG. 1: Comparison of the ground-state helium interatomic 
potentials obtained from the overlap theory, the extrapolation 
using the expansion of Ref.l24||. and ab initio calculations with 
exponentially correlated Gaussian functions [Tgj| 
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FIG. 2: Comparison of the ground-state neon interatomic 
potentials obtained from the overlap theory, the extrapolation 
using the expansion of Ref.|24|. with CI method [22^|. and in 
SCF calculations [23] 



ing good precision of both methods. The overlap calcula- 
tions agree with the expansion results over a larger range, 
with some deviation in the middle of the curve due to 
approximations of our theory. We can conclude that our 
method describes well the repulsive interaction in helium 
and the expansion method is suitable for calculations be- 
low 10 eV. In the case of neon (Fig[5J), we find close agree- 
ment between different methods: the overlap theory, 
the extrapolation from experiment |24| . configuration- 
interaction (CI) method [23], and self-consistent field 
(SCF) calculations [23| below 8 eV. We also found agree- 
ment between our method and that of Ref. |24| in argon, 
although not as close as in neon. Therefore, the ex- 
pansion of Ref. [24| gives the energies of ground state 
noble-gas dimers with sufficient precision for calculations 
of emission spectra and the analysis of self-absorption. 
The emission spectra of helium and neon molecules 



3 



1.0 - 




i i i i i i i i_ 

10 12 14 16 18 

Energy, eV 



FIG. 3: Theoretical helium and neon VUV emission spectra. 
The intensity is normalized to the maximum 

are shown in FigEI The initial He excited state has 
equilibrium distance R e = 1 A, energy 18.2 eV and the 
harmonic-oscillator frequency uj e = 1861 cm^ 1 17]. Cor- 
responding values for the Ne 3s 3 E„ initial state are: 1.75 
A, 16.11 eV, and 540 cm" 1 [Hill]. Neglecting anhar- 
monic corrections, we can write the nuclear wavefunction 
of the lowest vibrational state as 

a r a 2 (R-R e ) 2 ^ , . 

*( R ) = ^Tj2 ex P[ 2 1 (!) 

here a — (2nfj,cMj e /h) 1 ^ 2 , fi is the reduced mass, c is the 
speed of light. Although excited vibrational states are 
always populated in the process of any non-radiative ex- 
citation, very fast collisions and quenching will transfer 
much of their population into lowest vibrational states. 
Since only direct transitions are the most probable, 
the absorption spectrum will have the weighting factor 
^(i?)! 2 . The helium theoretical molecular spectrum is 
in better agreement with the emission spectrum observed 
in condensed helium discarge [3J, |25j than with that ob- 
served in liquid . High-pressure discharge emission 
has maximum located at 15.3 eV 0,113, while theoretical 
prediction is 15 eV, and the maximum of the emission of 
liquid is located at 15.5 eV and 15.8 eV 0. The width 
at half maximum of discharge spectrum is 2 eV 0, the 
same as our prediction, while the liquid spectrum has the 
width about 3 eV 0,0. From this comparison, we can 
conclude that the liquid effects shift the energy by 0.5-0.8 
eV and broaden the spectrum by 1 eV. A similar situa- 
tion occurs in neon. High-pressure neon discharge gives 
the maximum of emission at 15.5 eV, while the liquid has 
the maximum at 16 eV J]; the theoretical prediction is 
15.2 eV. The theoretical width is about 1 eV, and width 
of the liquid spectrum obtained with 1 eV resolution is 
about 2 eV. The liquid effects result in shift of at least 
0.5 eV and broadening less than 1 eV. The shift is not 
easy to explain; we can suggest that the interactions of 
more than two atoms can be responsible and the distance 
between an excited atom and closest ground state atoms 



can be different than that in the diatomic molecule. 

Apart from calculations of emission spectra, we also 
investigated various absorption mechanisms and calcu- 
lated self-absorption coefficients. In liquids as well as 
in dense gases in the vicinity of an atomic or molecular 
resonance collisional broadening is very important, and 
the absorption due to collisionally broadened Lorentzian 
wings is expected to be very strong; however, the impact 
approximation from which the Lorentzian shape-factor 
can be derived is not valid for far-off-resonance transi- 
tions, which is the case in our system, and instead more 
appropriate is to use quasi-static approximation (see for 
example |2^, pp. 241-245) in which the absorption pro- 
file is determined from probabilities W(R)dR of nearest 
particles within the range of distances (R, R+dR) from a 
given atom. Either the classical distribution or quantum- 
mechanical wavefunctions can be used to obtain these 
probabilities. In the classical case, the Boltzman dis- 
tribution for low temperatures of the liquids gives very 
strong exponential suppression for the detuning of or- 
der 1 eV and the absorption can be neglected. On the 
other hand, the consideration of wavefunctions of nuclear 
motion (the tunnelling into classically-forbidden region) 
gives negligible increase in absorption because the tun- 
nelling probability is small and the energy deficit will re- 
main whether the tunnelling occurs in the ground state, 
excited state, or both. Finally, we conclude that the dom- 
inant absorption is likely due to wings of the radiatively- 
broadened profile. 

For large detuning and radiative broadenin g, t he ab- 
sorption coefficient a a is (see for example Ref.[27|) 

_ Ne 2 uj jojf 

aa ~ II 9\? i 2 2 \ I 

cmeo (oJq — u> z ) A + ^ z io z 

where 7 is the width of the resonance, ujq is the resonance 
frequency, uj is the frequency of the light, N is the num- 
ber density, c is the speed of light, en is permittivity con- 
stant, / is the oscillator strength of the dominant atomic 
transition. The equation assumes CI units. Note that 
the frequency dependence approaches Lorentzian for de- 
tuning much smaller than the resonance frequency. The 
self- absorption coefficients for liquid helium and neon cal- 
culated using the above equation are given in Table [I] 

The scattering lengths can be found from the theory 
of Rayleigh scattering given by Landau and Lifshitz [28| : 

We neglected the temperature derivative term which was 
estimated to be small in Ref.0. The compressibility Kt 
of He and Ne are 208xl0~ 10 and 4.95xl0~ 10 cm 2 /dyne; 
p is the density. The dielectric constant is related to the 
refractive index: e = n 2 . The refractive index can be 
calculated using the equation (see for example Ref. |2jJ) 

n 2 -l = a 2 Na X 2 v A 2 /, 

n 2 + 2 3tt ^ X 2 - A 2 U 
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TABLE I: Liquid helium and neon absorption and scattering 
coefficients. Parenthesis denote of the uncertainty of calcula- 
tions in the last digit. 



Liquid 


n, cm 3 


\,nm Xo,nm 


a a , m" 


-1 


a s , m 1 


He 


1.88 xlO 22 


78 58.43 


0.36 


0.028(2) 


0.21(4) 


Ne 


3.60 xlO 22 


80 73.59 


2.02 


0.15(4) 


1.5(8) 



where a is the fine-structure constant, N is the number- 
density of atoms, and ao is the Bohr radius. The summa- 
tion is extended over continuum. For neon, we included 
9 values from Ref . |30| and estimated remaining contri- 
bution from the Thomas-Reichc-Kuhn rule. For helium 
we calculated 40 oscillator strengths with CI method us- 
ing B-spline basis sets obtained in a cavity 30 a.u. The 
continuum contribution in helium is larger than that in 
neon because of larger detuning from the dominant tran- 
sition for the wavelengths of interest. To check our con- 
tinuum contributions in helium, we calculated the refrac- 
tion index at wavelengths in the range from 168 to 288 
nm to compare with available the precise measurement 
by Smith et al. j2|j. We found that our ab initio val- 
ues agree well with the experimental values. The largest 
uncertainty in our helium calculations at 78 nm comes 
from the replacement of continuum states with quasicon- 
tinuum states, and in neon calculations from the estima- 
tion of continuum contributions using the sum rule. The 
results of calculations of scattering coefficients are given 
in Table IJ 

We found that our results for scattering lengths and 
dielectric constants, e# e = 1.086±0.006 and ejv e = 1.5± 
0.2, are in agreement with those obtained by Seidel et al. 
[{| using extrapolation of experimental refraction indices: 
(He = 1.077, ejv e = 1.52, Rh £ — 600 cm , and Rnc = 60 
cm. Our calculations has been done to check the previous 
results and to evaluate the uncertainty. 

In our calculations we found that absorption lengths 
are smaller than the size of the neutrino detector, so 



that it is important to know whether the excited atom 
re-emits an absorbed photon or it is quenched non- 
radiatively. Because for non-radiative decay triplet and 
singlet states are similar and in Ref.|3lj| it has been found 
that lifetime of the helium triplet state in solid exceeds 
15 s, we can conclude that quenching is a slow process 
and the total emission power will not be reduced much. 
The absorption-re-emission process is equivalent to the 
scattering in this case. 

In this paper we have analyzed UV emission, self- 
absorption, and scattering of liquid helium and neon. 
For emission spectra we have found reliable data for he- 
lium and neon molecules and developed our own model 
of repulsive interaction to check and understand the ex- 
isting calculations. Using potential curves we obtained 
emission spectra which agree well with experimentally 
observed condensed discharge spectra, although liquid 
emission is shifted due to liquid effects. We have con- 
sidered various mechanisms of self-absorption and found 
that the radiative broadening is the dominant mecha- 
nism if the emission absorption offset exceeds 1 eV. The 
self- absorption is significant at the scale of the proposed 
neutrino detector and can impose restrictions on the ex- 
periment; fortunately, non-radiative relaxation seems to 
be slow and the absorption results mostly in re-emission 
so that it is similar to scattering in terms of conservation 
of the number of photons, but dissimilar to it in terms of 
photons' propagation dynamics. We have also calculated 
dielectric constants and Rayleigh scattering parameters. 
Results are in agreement with previous calculations. We 
understand now dominant effects of the emission, self- 
absorption, and scattering, but some questions remain 
and further investigation is necessary: the study of the 
liquid shifts and other possible liquid effects, the study 
of quenching mechanisms, the measurement of scattering 
and absorption lengths. The author is grateful to Dr. 
McKinsey for valuable discussions and comments on the 
manuscript. 
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